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The quantum electrodynamic theory of the nuclear recoil effect on the atomic g factor to all 
orders in aZ and to first order in m/M is formulated. The complete aZ-dependence formula for 
the recoil correction to the bound-electron g factor in a hydrogenlike atom is derived. This formula 
' is used to calculate the recoil correction to the bound-electron g factor in the order (aZ)'^m/ M for 

' an arbitrary state of a hydrogenlike atom. 

(N ■ 

^ ' I. INTRODUCTION 

High precision direct measurements of the bound-electron g factor in \aw-Z hydrogenhke ions, which are performed 

(~| ■ by a GSI - Universitat Mainz collaboration , have triggered a great interest to the calculations of QED and nuclear 

CIh] corrections to this effect . It is expected that in nearest future the collaboration will extend these measurements 

_L( to high-Z systems. The accurate measurements of the g factor in heavy ions would provide a good possibility for a 

G test of the magnetic sector of QED in a strong Coulomb field. To date, investigations of the bound-electron g factor in 

O , heavy ions were performed only by an indirect method which is based on using the fact that the bound-electron 

d ' g factor can be expressed in terms of the transition probability between the hyperfine splitting components Jl^ , p^ . 

^ , In particular, in the case of hydrogenlike bismuth, it was found that including the QED correction to the g factor is 

O ■ necessary to obtain agreement between theory and experiment (see [ p"5|jl^ ] for details). However, the indirect method 

' ^ ! '^ot provide so high accuracy as the direct method if the latter is extended to high-Z systems. 
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It is known (see, e.g., 17|) that calculations of heavy ions have to be performed without any expansion in the 
' papameter aZ. In particular, it means that, to zeroth-order approximation, the electron in a hydrogenlike atom must 
be described by the Dirac equation with the Coulomb potential induced by the nucleus. For the point-charge nucleus. 



the bound-electron g factor in a hydrogenlike atom is easily calculated analytically (see, e.g., 18|). The one-loop 
QED corrections to the Is g factor were calculated to all orders in aZ in (see also 0). The finite nuclear size 

^ correction can be found by solving the Dirac equation for an extended nucleus For low- and middlc-Z atoms, 

■ this correction can also be found by non-relativistic and approximate relativistic formulas 191. The nuclear recoil 

o 



correction to the Is g factor was evaluated only to the lowest order in aZ [p0|-p3t (see also [ ^4| , p^ ). The methods 
developed in [po|-p3| are suitable for calculations to a few lowest orders in the parameter aZ and are not intended 
, for calculations of the recoil effect to all orders in aZ. Therefore, first of all we need to develop a new method which 
would provide a systematic treatment of the recoil effect on the g factor to all orders in aZ and to first order in 
, m/M (m is the electron mass and M is the nucleus mass). Such a method is proposed in the present paper. This 
1 method is valid for the general case of a many-electron atom. ApUieng the method to a hydrogenlike atom, the 
O complete aZ-dependence formula for the recoil correction to the g factor is derived. To the lowest relativistic order 
'U^ ('^^ (ctZym/M), this formula is used to calculate analytically the recoil correction to the g factor in the case of an 
arbitrary state of a hydrogenlike atom. For the Is state, the result obtained in this paper coincides with the related 
: result obtained previously in [p0|-p3[. 
O-i The relativistic units (?i = c = 1) and the Heaviside charge unit {a = e^/(47r), e < 0) are used in the paper. 



II. QUANTUM ELECTRODYNAMICS TO FIRST ORDER IN m/M BEYOND THE FURRY 

PICTURE 

To derive the Hamiltonian of the atom to first order in m/M and to lowest order in a (but to all orders in aZ), 
we start with the relativistic quantum mechanical treatment of the electrons and the non-relativistic treatment of the 
nucleus. For simplicity, we consider the nucleus as a spinless particle. The electrons and the nucleus interact with each 
other, with the quantized electromagnetic field, and with the classical homogeneous magnetic field, Aci(r) = [Hx r]/2. 
In the Coulomb gauge and the Schrodinger representation, the Hamiltonian of the system to first order in the 
interaction with Aci(r) is 

N' 

H = ^[a, • p(^) + ftm + V{r['^ - r(")) - ea, • A(r(^))] 

1=1 
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-ef:«..A.,(rr))-^p(").Ae,(r(")), (1) 

i=l 

where N' is the total number of the positive and negative energy state electrons, p^-'^'' and p^"^ are the electron and 

nucleus momentum operators, respectively, V^(r|'^'' — r*^"^) is the operator of the electron-nucleus Coulomb interaction, 
A(r) is the vector potential of the quantized electromagnetic field, 7i = V x A, and £t = —dA/dt. The last two 
terms in equation (|l|) describe the interaction of the electrons and the nucleus with the homogeneous magnetic field. 
We introduce the center-of-mass variables for the electron-nucleus subsystem 

1 ^' 

R= fMr(")+myrl'^V (2) 

M + N'm V ^ y ' ^ ' 

fc=i 

r, = > - r(") . (3) 



From these equations we derive 



rr' = + R Tfc , (4) 

* A/I A- N'm /-^ ' ^ > 



N' 

M + N'm 

fe=i 

N' 

r(")=R !!1 yr„. (5) 

k=\ 



For the corresponding momentum operators (P = — ^Vr, Pi — iVrJ, we obtain 

(e) _ m 



M -I- N'ra ^-^ 



k=l 



To keep in the Hamiltonian the terms of zeroth and first orders in m/M, we can replace equations (^)-(|^) by the 
following ones 



N' 



M 

fe=l 





= R 


N' 

m \, ^ 

' M ^ 

fc=i 


(e) 
PI 


= P. ^ 


mp 

M ' 

N' 


pin) 


= p 


k=l 



(9) 
(10) 

(11) 

With these substitutions, the individual terms in equation (|^) are transformed as 

N' N' N' 

^a,.pfU^a,.p, + ^^a,.P, (12) 

2 — 1 2—1 i—1 

N' N' 

^y(r(^)_rH)^^y(r,,), (13) 



i=l i=l 
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N' N' N' 

-e^ a, • A(rf = -e^ a, • A(r, + R - ^ 



i=l 



1=1 
N' 



k=l 



-e^ai- A(r, +R) 



1=1 
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k=l 

a 



2 ^ l'"^ - 



-S[(i:-^)(i:''.-A(o 

IE 



(14) 
(15) 



— [pH ~ |e|ZA(r("))]^ = ^ - • E + ^ (E P^- 



fc=i 



2Af 



fc=i 



— P • A R + 



N' 



k=l 



5:p,).a(r) + ^a^(r) 



(16) 



AT' N' 

-e^OLi- Aci(r|''^) = -e^ CKi • Aci(^ 



r, + R 



i=l 



M 



fe=l 



JV' N' 

-e^ai- Aci(ri) - e ^ • Aci(R) 



1=1 1=1 

TV' Af' 

+^^E"^ ■E^^i(^^)' 

i=i fc=i 



^p(")-Aei(r(")) 



JV' 



^I^p.a.kR) + ^Ep^-AcKR). 



Af 



(17) 
(18) 



fc=i 



Here we have disregarded the terms of second and higher orders in m/M. Because the second term in the right hand 
side of equation (|lj) contains the factor e-^, to first order in m/M, it can contribute only to first and higher orders 
in a. Therefore, to the order under consideration (to zeroth order in a, to first order in m/M, and to all orders in 
aZ) this term can be omitted. We obtain 



AT' ^ 

H = ^[a, ■ p, + Am + V(v,)] + t^YI 



2— |r.-r.| 

N' 



+ 1 [ [^'(x) + W'(x)] - e a, • A(r, + R) 

1=1 

+ 2M + 2MlEP'^j +mE"-P 



1 ^' 

P Vpfe 

M ^ ^ 



/c=l 
JV' 



A/ 



P A(R) 



-^(i:p^)-a(r)+^:^a^(h) 



k=\ 



JV' JV' 

-e ^ • Aci(ri) - e ^ a,; • Aci(R) 



JV' JV' II 
™ \ ^ \ - . , , \&\Z ^ 
+^ME"^-EAci(rfc) 



i=l fc=l 



Af 



P-Aei(R) + ^5]p;.-Aci(R) 



JV' 



(19) 



fc=i 



3 



In the second quantized form, this Hamiltonian can be written as 

13 

s=l 



H^Ho + J^H^l (20) 



Ho= Jdx V'^(x) HV ■ a + /3m + V (x)]V^(x) ~ ^ J Ax0(x) 

+ 1 J dx[f2(x)+7i2(x)], (21) 



H^nl = -e y dx dy (x)aV(x) A(x + y)0t (y)0(y) , (23) 
= ^ / '^'^ (x) (-^ V.) V^(x) 1 dy (y) (-» Vy ) V(y) , (24) 
^iift^ = [d^ V^t(^)(_,Vx)^(x) / dy A(y)0t(y)0(y) , (25) 



H^S -^Jdy A'{y)^Hymy) , (26) 

H^n! = ^ /^^ ^t(x)«V(x) 1 dy 0t(y)(_,Vy)0(y) , (27) 

H^I! Vt(x)(_jVx)^(x) 1 dy </)t(y)(-zVy)0(y) , (28) 

= / A(y)0t(y)(_,Vy)<^(y) , (29) 



rr(9) 
^int 



-e j dx?At(x)Q,. Aci(x)V'(x), (30) 



^iLf = 2M / '^^ V'^(^)"V'(X) / dy ^t(y)A,,(y)^(y) 



+~ j V^^(y)Aci(y)^'(y) y dx v/(x)aV'(x) , (3i) 

^-e I V^t(x)a^(x) / dy 0t(y)A,i(y)0(y) , (32) 



^il?^ = / rfy A,i(y)0t(y)(_^Vy)0(y) , (33) 



= ^ / '^xV^t(x)(_,Vx)^(x) j dy 0t(y)A,i(y)</>(y). 



(34) 



Here V'(x) and (/>(x) are the electron and nucleus field operators in the Schrodinger representation and Ax is the 
Laplace operator. The scalar product is imphcit in equations (p3|)-(p^. It should be stressed that, despite the nuclear 
field operator 0(x) is introduced, the Hamiltonian (20) has physical sense only in the one nuclear particle sector. The 
introduction of 0(x) will simplify the use of the standard QED methods for calculations by perturbation theory. 
The total momentum operator of the system is given by 

P = y"dx (/.t(-iVx)0(x) + j [ft(x) X 7^(x)] . (35) 

This operator commutes with the Hamiltonian of the system if the external magnetic field is switched off, Aci(x) — 0. 
It also commutes with the operator 

s=W 

H' = H,+ Y^ H^:^ . (36) 

s = l 

In the theory without the terms H^^^P , H^nt \ ^'^'^ ^inf'' ^ ^ i'^ the theory without the external magnetic field 
at all, operators (/)(x), A(x), and '(/'(x) obey the following space-shift transformation rules 
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(/)(x + a) = exp (— ia • V)(j){yi) exp (ia • V) , (37) 

A(x + a) = exp (-ia • P)A(x) exp (ia • V) , (38) 

'0(x) = exp (— ia • 7^)-0(x) exp (ia • V) . (39) 

For any operator -F'(x), the transition to the Heisenberg representation is performed by the substitution 

Fs{x) = exp {-iHt)FH{t, x) exp {iHt) . (40) 



The calculations of the atomic energy levels, based on the Hamiltonian (^0|), can be performed by perturbation 
theory. It is convenient to consider that in zeroth approximation the Hamiltonian of the system is given by the term 



Hq. The terms H^^^ - H^^^ are accounted for by perturbation theory. As usual (see, e.g., [^5y26[ ), the vacuum state is 



defined as a state in which the negative-energy continuum is occupied by electrons and the electron current operator 
eip{x)^^ip{x) must be replaced by (e/2)['0(a;)7^, ^(x)] in all the equations given above. 

III. GENERAL FORMULAS FOR THE ENERGY SHIFTS 

First, we formulate a procedure for deriving the energy shift due to the terms - h[^\ The operator H' defined 
by equation (^) commutes with the total momentum of the system ( p5| ) and, therefore, the atom can be characterized 
by a fixed momentum K. We are interested in the energy shift, Ai?(a k) = ^(a.K) ~ -^(a^K)' "-"^ ^ single isolated level 
|a, K) of an TV-electron atom. Here N denotes the number of the atomic (positive-energy-state) electrons, K is the 
total momentum of the system described by the Hamiltonian H' , and a denotes the set of the other quantum numbers. 
In zeroth approximation, the energy £'1^''^) 

^^(alK)=£a,+---ea„+^. (41) 

The one-electron energies e„ are defined by the Dirac equation 

(-ia • V + /3m + V^(x))V'„(x) = e„V«(x) , (42) 

where l^(x) is the Coulomb potential of the nucleus. We note that the formalism considered here allows us to 
account partially for the finite nuclear size effect by employing in (|4^) the potential of an extended nucleus. We 
could start also with the Dirac equation with an effective potential V|,ff(x) (e.g., a local version of the Hartree- 
Fock potential) which approximately describes the electron-electron interaction. In this case the interaction with 
the potential Ay(x) = l^(x) — T4ff(x) must be accounted for perturbatively to eliminate double accounting the 
interclcctronic interaction corrections. The electronic part of the unperturbed wave function is a linear combination 
of one-determinant functions 

ua(xi, XA,) = ^ C^-L. Y^^^i)P^p^^ (xi) . . . ^p,„ (x^,) . (43) 

b ^ p 
We introduce the Green function ga{E,'K.) by 

g,{E,K.)5{E' ~ E)5{K.' -K) ^ ^^-^ [ dxi • • • dxArdx'i • • • dx^dXdX' 

/oc 
dt dt' exp {iE't' - iEt) exp (-iK' • X' + iK • X) 
-OO 

xwt(x;,...,x'^)(o|rV'(i',x;)...V'(t',x'^)0(i',x') 

X(/)t(t,X)?At(i,xjv) • • • V^(t,xi)|0)ua(xi, , (44) 

where V'(^,x) and 0(t, x) are the electron and nucleus field operators in the Heisenberg representation and T is the 
time-ordering operator. As one can see from equation (^), ga{E, K) is the Fourier transform of the two-time Green 
function. The ga{E,Vi) function contains the complete information about the energy levels of the system described 
by H' . Defined by equation (jij) for real E, it can be continued analytically to the complex E plane with some cuts 
(see, e.g., [ p7|p^ ). From the spectral representation of ga(£',K), we obtain 

gaiE, K) = ^ — — - + terms that are regular at £' E^^ k) i (45) 

E — E(^a,K.) 
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where 

/ dxi...dxAr«t(xi,...,xjv)(O|V(O,Xi)...V(O,XA,)0(O,O)|a,K 



ip^a.K) = ^^=- J dxi---dxAr«I(xi,...,xjv)(O|V'(O,Xi)---V(O,XA,)0(O,O)|a,K>. (46) 
As in , using equation (^5|) , we can derive the foUowing formula for the energy shift 

1 /".„,„ -^(0) 



Ai^(a,K) = , . , (47) 

l + —idE Ag,(£;,K) 

where l^ga{E,K) = ga{E,K) - g^a\E,K) and g^a\E,K) ^ {E - ^he contour T surrounds the pole 

corresponding to the level a and keeps outside all other singularities of ^^(i?, K). It is assumed that the contour T 
is oriented anticklockwise. The Green function ga{E^'K.) is constructed by perturbation theory after the transition in 
equation ( ^4|) to the interaction representation. Since we are interested in the energy shifts in the atom rest system, 
we must put K in formula (|47|). 

Formula ( ^ ) can be used to calculate the corrections due to the terms hI^-HI^^\ The term Ai/', defined as 

Ai/' ^ Hl^ + + , (48) 

does not commute with V and, therefore, requires a special treatment. To first order in the interaction with the 
magnetic field, the energy shift due to this term can be written as (cf., pof) 

Ai?(,,K)'5(K' - K) ~ (a, K'\AH'\a, K) , (49) 

where |a,K) and |a,K') are the eigenvectors of the operator H' . Since we are interested in the corrections of first 

order in the magnetic field, in equation ( ^ ) we can consider the eigenvectors of the operator H" = H' — H^^^ — H^^^^^ . 
We assume that they are normilized by 

(a,K'|a,K> =(5(K'-K). (50) 
We introduce the Green function ga{E' , E, K', K) by 

ga{E', E, K', K) = Jd^i--- d^Nd^[ ■ ■ ■ dx'^dXdX' 

/OO 
dt dt' exp {iE't' - iEt) cxp {-iK' • X' + iK • X) 
-OO 

xut (x;, x^)(0|TV'(t', x'l) • • • ^(f, ^Mt', X')Ai/'(0) 

x0t(t,X)V^t(i^xjv) • • •^^(t,Xi)|0)Ma(xi, ...,xjv) , (51) 

where, as in equation (^), the Heisenberg representation is used and Aff'(O) = AH'{t = 0). The analytical properties 
of this type Green function as a function of two complex variables E and E' in the region E ~ £'(a,K)j E!' E(^g_^j^r 



can be investigated by deriving the double spectral representation (see, e.g., |28|j29|]). From this representation, we 
obtain 

g,{E\E,K\K) = (a,K'|Ag-(0)|a,K) /y) 

+ terms that are regular functions of E' or E when 

E' - £;(a,K') and£; ^ E^a,K) ■ (52) 
As in the case of the photon emission by an atom considered in detail in p8| , from equations ( |49| ) , ( |5^ we derive 

J_IdE^(fdE' ga{E', E, K', K) 

A£;(.,K)^(K' - K) ^ f ^ ' (53) 

— j^^dE ga{E,K')\ [— j>^dE ga{E,K) 

where the contours F and F' surround the poles corresponding to the levels £'(a,K) and -E'(a,K')j respectively, and 
keep outside all other singularities. The Green function ga{E' , EjK.' ,K.) is constructed by perturbation theory after 
the transition in ( pl| ) to the interaction representation. Since we are interested in the energy shifts in the atom rest 
system, at the end of the calculation by formula (El), we must put K = 0. 
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IV. RELATIVISTIC FORMULA FOR THE RECOIL CORRECTION TO THE 
BOUND-ELECTRON g FACTOR IN A HYDROGENLIKE ATOM 



The method formulated in the previous sections can be used to derive the complete aZ-dependence expressions for 
the recoil corrections to the energy levels (Aci(x) — 0) and to the atomic g factor to first order in m/M. The relevant 
expression for the recoil correction to the energy levels in a hydrogenlike atom was first derived by a quasipotential 
method i n po| and was subsequently rederived by other methods in (a part of the expression was previously 

found in |]34|). In this section, we derive the corresponding formula for the recoil correction to the bound-electron g 
factor in a hydrogenlike atom. 

In what follows, we will consider that V{x) in the Dirac equation (|4^) is the potential of an extended-charge nucleus. 
This will allow us to avoid a special treatment of the degenerate states of different parity (such as 2s and 2pi/2) in the 
derivation of the contributions from the AH' term. However, as it will be seen from the final formulas obtained below, 
they have neither singularities nor ambiguities when ^(x) — > Vc(x) — —aZ/\x.\ and, therefore, the pure Coulomb 
potential can be used in these formulas as well. 

We are interested in the energy shifts of first order in the magnetic field and of first order in m/M. The contributions 
of interest arise in a few lowest orders of the perturbation theory. They can be conventionally represented by the 
following combinations: 

+hI!I X Hl^i X hII^ X Hill + permutations , (54) 
A£;(^")~i?(f , (55) 

AE^'"' - Hl^' X His + Hi:^' - Hl^ 

+hIII'> X hI^I X hI^I + permutations . (56) 

(12) (13) 

As to the terms H^^^ and H^^^ , in the order under consideration, the first of these terms gives a contribution which 
is independent of the internal atomic quantum numbers while the second one gives zero when it is averaged with the 
electron wave function of a definite parity. Therefore, the both terms can be omitted. 

To calculate AE^^\ it is convenient to consider that H^^ is included in the unperturbed Hamiltonian. This means 
that the interaction with the magnetic field, (5V^(x) = — ea • Aci(x), must be included in the Dirac equation (|4^). 
Then the derivation of the formula for this correction reduces to the related derivation of the recoil correction to the 
energy levels for Aci(x) = |^-|||. It follows that to obtain AE^'^\ we may simply replace V^(x) by V^(x) -I- 6V{x.) 
everywhere in the formula for the recoil correction to the energy level and decompose it to first order in 6V{'x.). 
However, before to do that, let us demonstrate how the recoil correction to the energy level can be derived within the 
formalism considered here. 

Let us consider, for example, the contribution ~ ^^i'nt x + H^ x ij/^j . In the interaction representation, the 
related contribution to the two-time Green function is 

AG'(t',x',X';i,x,X) = ^-^ / dtidt2 (0|TV'(t',x')</'(t',X')V'Ut,x)(/)t(t,X) 

x{^^ / '^''I'^^i [V^(ii,xi)a\^(ti,xi)]v4^(ti,xi+Xi) 
x0t(t,,Xi)0(ii,Xi)^ y dyidYi ^Ht2,yi) 

X (-zV^ J^(t2, yi)A'=(t2, Yi)0t (t2, Yi)0(t2, Yi) 
\e\Z 



M 



J J dyidYi V't(i2^y^)(_,v^j^(i2,yi)A^-(^2,Yi) 
xcj)^{t2,Y,)cl>{t2,Y,)^ J dxidXi [V't(ii,xi)a\^(ti,xi)] 
xA^(<i,xi+Xi)0t(ti,Xi)0(ti,Xi)}|O)co„, (57) 

where all the operators are considered in the interaction representation and the summation over the repeated indices 
(i, A: = 1, 2, 3), which enumerate components of the three-dimensional vectors, is implicit. The label "con" means that 
contributions containing disconnected vacuum-vacuum terms must be omitted. According to the permutation rules 
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for the T product of the boson and fermion operators, second term in the right hand side of equatio n (]57| ) is equal to 
first one. Assuming that the time-ordering operator is defined for equal-time fermion operators by ||l7| 



T[A{t)B{t)] ^ \A{t)B{t) \B{t)A(t) , (58) 



we have 



I 1 7 

AG(i',x',X';i,x,X) = / dtidiz (O|T?/.(t',x')</>(t',X')?/'t(t,x)0^(t,X) 

X j dxidXi V^(ti,xi)aV(ii,Xi)A'(ti,xi +Xi) 
x0t(t^,Xi)0(ti,Xi) j dyidYi Vt(t2,yi) 

X {-iVlJij{t2.y^)A\t2, Yi)0t (^2, Yi)0(t2, Yi)|0)co„ . (59) 

Using the Wick theorem and keeping only the terms which will contribute to the state-dependent energy shift, we 
obtain 

AG{t',x',X';t,x,X) f dtidt2 f dxidXidyidYi[S{t' -ti,x',xi) 

X a'S{t^ - <2 , xi , yi) (-iV^ J5(t2 - yi , x) 
+S(t' - t2, x', yi)(-iV^J5(t2 - ti, yi, xi) 
X a'Sitr - t, xi,^)]D'^ ih - ^2, xi -f Xi - Yi) 
X [Snif - ti,X' - Xi)5„(ii - i2, Xi - Yi) 
xS„{t2 - i, Yi - X) + Sn{t' - t2,X' - Yi) 

xS„{t2 -h,Yi- Xi)5„(ti - i,Xi - X)] . (60) 

Here 

GW + ^ * ^ r ■ ^mV^ V'n(x)V4(y) 

S{t-t,x,y)^— duj exp[-iuj{t ~t)]2_^ — — 61) 

is the electron propagator in the nuclear potential V^(x), 

is the transverse part of the photon propagator in the Coulomb gauge, and 

, i , , , f dK exp [iK • (X - Y)l , , 

Sn{t'-t,X-Y) = — duj cxp [-iuj{t' - t)] / -—3 (63) 



27r7_"- '"J (2^)3 c.-lii+zO 



is the nucleus propagator. In equation (|6^) the index n runs over all bound and continuum states. Since we are 
interested in the corrections of first order in m/M and we have already the factor 1/M in front of expression (|60|), we 
can consider the limit il/ — s- in the expression for the nucleus propagator. In this limit, we have 

Sn{t' -t,X-Y)=e{t' -t)6{X-Y). (64) 

Using this expression for S'„, we obtain 

AG(t',x',X';t,x,X) = j dhdt2 J dxidyi [S{t' - h,x' 

xa'Sih ~ i2,xi,yi)(-iVjj5(t2 - i,yi,x) 
+S{t' - t2, x', yi)(-iV^ J5(<2 - ti, yi, xi) 

XQ!*5(tl - t,Xi,x)]L»*''(ti - t2,Xl) 

x,5(X' - X)[9{t' - h)9{h - t2)0{t2 - t) 

+e{t' - t2)9it2 - h)9{h - t)] . (65) 
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The related contribution to the Green function (7a(£',K) defined by equation ( }44D is 

Aga(E,K)S(E' - E)5(K' -K) = -^-^^^ f dtdt'dtidt2 / dx dx' dX dX' 

27^^ [2ny M J 



xdxidyi exp {iE't' — iEt) 

X exp {-iK' -X' + iK- X)ijl{x')[S(t' - ti, x', xi) 
xa*S'(ti -i2,Xi,yi)(-iV^JS'(<2 -i,yi,x) 
+Sit' - h, x', yi)(-zV^ JS-lia - ti, yi, xi) 
xa'^(ti -i,xi,x)]?/;a(x)i?'^(ti -i2,xi) 
x,5(X' - X)[0(t' - h)e{h - t2)0(i2 - 1) 

+e{t' - t2)9{t2 - h)e{ti ^ t)] . (66) 



Integrating over X, X', and over the time variables and setting K = 0, we obtain 

e\e\Z 1 Z""" , fv- {a\a'-D'^{u;)\n){n\{-iV'')\a) 



^ M [E^SafJ^^ £;_^_e„(i_^o) 

-iV'^)|n)(n|a^i:»'''''(c 
£;-cj-e„(l-iO) 



(a|(-»Vfe)|n)(n|a'jj"--H|a) ^ 
^ i^-c.-e.„fl-zO) 1' ^''^ 



where 



^ (2^ .^ k^+.o (^"-l^) 

= _^|£^P(iW^^^ ^ V'V'(°"P(^y-^H . (68) 
47r L r w^r J 

The corresponding contribution to the energy shift according to equation (^7|) for K = is 

= 1^ A r dojS^S {a\a^D^H^)\n){n\{-iV'')\a) 

M 977 / 



M 27ry_^ ^[ Ea - LV - en{l - iO) 

{a\{~iV'')\n){n\aW''{uj)\a) 1 
- iO) )■ 



(69) 



The energy shifts due to the terms H^^l and H^^^ (~ H^^l x H^^l x i/j^jj^ + permutations) can be derived in the same 
way. The total contribution of the terms H^^ , h[^1 (^ h[^1 x H^H + x h[^1), and TJ^f,' (^ TJ^f,' x TJ^f.^ x H^^ 



permutations) is 



where p'^ = — iV* 



and 



AE=^ I dto (a|b^ - i?'^(a>)]G(c. + £,)[/ - i?^(c^)]|a) , (70) 



D^iuj) = -AnaZa'D^^iuj) , (71) 



^c^-e„(l-zO) 

Deriving equation ( |70| ) we considered that in the zeroth approximation the electron obeys the Dirac equation with 
the nuclear potential V{x). As was indicated above, to obtain the desired formula for AE'^^\ defined by equation 
(P^), we should replace F(x) by V{x.) + SV{x.), where 6V{x) — —ea ■ Aci(x). in equation ( [70| ) and expand it to first 
order in SV{'X.). As a result of this expansion, we obtain 
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+(a|[p'= - D''{u)]G{lo + e„)[p'= - D''{u)]\5a) 
+(a|[/ - D^{uj)]G{uj + ea){5V - Ssa) 

xG{uj + ea)[p'' - D\io)]\a)} , (73) 



where 



5ea = (a|JT/|a) , (74) 
|(5a) = 2^ _ ■ (75) 



For practical calculations, it is convenient to represent expression ([73|) by the sum of the lower-order term and the 
higher-order term, 

Ai;(i) = A£;[^^ + A£;^^^ , (76) 
a4'' - ^{('^^It^'' - '2D\Q)]p^\a) + (a|[/ - 2D\Q)]p^\5a)} , (77) 



+ (a|(^^(.)-f^)G(. + e.)(..^-(.) + 5^)l^« 

+ - ^-^^)g{lO + ea)(<5F - <5ea)G(L^ + Sa) 

\ UJ + l\J / 

x(i-'H+^''''l' 



-H(a|(D^.)-mG(. + e.)i^|a)}, (78) 

where [A, B] = AS - BA. The term A£'{^^^ contains the non-relativistic and lowest-order relativistic contributions 
and a part of the higher-order relativistic contributions. The term ^e'^^ contains all the higher-order relativistic 
contributions which are not included into AE'^^''. 
The direct evaluation of the term hI^^^ yields 

du; [{a\a^G{u;)A';,\a) + {a\Al,G{^)a^\o)] ■ (79) 



The calculation of the correction AE''^''^ is a more difficult task. Here we derive in detail the contribution which 
appears as a combination ~ -ff/nt^'' x ^^i'nt ^ -f^/nt + permutations. To find this correction we should use formula ([s^). 
The related contribution to the Green function ga{E' , E, K', K) defined by equation ([5l| ) is 

1 2 1 I 7 /'OO r 

Aga{E\E,K\K) - — dtdt' dtidt2 /rfxdx'rfXdX' 
(27r) M J -oo J 

X exp [iE't' - iEt) exp {-iK' -X' + iK- X)7/>t (x') 
x(O|TV'(i',x')0(t',X') J <izidZi^^(0,zi)aV(0,zi)A^i(Zi) 

X(/.t(O,Zi)</.(O,Zi)0t(i,X)V/(i,x) J dxi rfXii^t(i^^xi) 
xa*i^(ti,Xi)A'(ii,xi -f Xi)0t(ti,Xi)(^(ti,Xi) 

X y"dYi ^'=(t2,Yi)0t(t2,Yi)(-iV^J0(i2,Yi)|O)eon^a(x), (80) 
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where we have used equation (p8|). To evaluate this expression we employ the Wick theorem and equation for 
the nucleus propagator. Keeping only the terms which will contribute to the state-dependent energy shift, we obtain 

Aga{E\E,K',K) = -^-^^ J dtdt'dtidt2 J dxdx'dxirfzi exp {iE't' ^ iEt) 

X [i^l{K')S(t' - <i, x', xi)a^5(ti, xi, zi)a'5(-t, zi, x)^,(x) 
+Vl(x')^(i', x', zi)a'5(-ti, zi, xi)a^5(ti - t, xi, x)Vj,(x)] 
xD'\h - i2,xi)[0(i' - h)eih - h)e{t2)9{-t)B'^'{K',K) 

+e{t' - h)0{h)e{-t2)e{h - t)B^^{K\K) 

+e{t' ~ t2)9{t2 - tx)9{tx)e{--t)B'f (K', K) 

+e{t' - t2)e{t2)e{~h)0{ti - (k',k) 

+e{t')9{~t2)e{t2 - h)e{ti - t)B^\K',K) , (8i) 



where 



(K', K) = ^ J dX exp i-zK' ■ X)(-zV^)A^i(X) exp (zK • X) , 
B^-'(K',K) = ^ J dX exp(-zK'.XM^i(X)HV^)exp(zK.X). 



Integrating over the time variables, we obtain 

„2 



Aga{E',E,K',K) = -^^^ duj J^{a\aW^'' {oj)\n) {n\a' \a) 

x(0(£„) 1 1 1 —B1'{K\K) 



xB^''(K',K) + 0(£ 



E' -£a E' - U} - En + iO E - UJ - Ea + lO E - Ea 

i i 



E'-En E~£. 

i i 



E' -Ea E' + OJ - Ea + iO 

(K',K) 



E' -Ea E' +UJ -Ea+iQ E + UJ + En - 2Ea + lO E - Ea 



(K',K) 



E' -Ea E - UJ - Ea+M 



X B^'(K',K) 

E + En — 2£a E — Ea 

i i 



E' -Ea E + UJ + En ~ 2Ea + iO E + En - 2Ea E - Ea 



(K',K)} 



E' -Ea E' -UJ + En - 2Ea + iO E' + En - 2Ea E - Ea 



xi3f(K',K)-0(-e„)- 



E' -Ea E - UJ - Ea + iO 

—b!^'(K',K) 

E' -UJ + En-2Ea+iO E-Ea ^ ^ ' 



E' -Ea E' +UJ-Ea + iO E' + En - 2Ea E - Ea 



(82) 

(83) 
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xSf (K',K) + 6i(e ^ * 



E' -£a E' +uj-ea + iO 
(K',K) 



E + UJ — En + E — Ea 

„, s i i i i 



E' - Ea E ~ UJ - Ea + iO E - En E - Ea 

i i 



xB^''(K',K) + 6i(e, 



i?2 



E' ~Ea E + LJ - En + iO 

'='(K',K)}. (84) 



Taking into account that (n|a'|a) = if e„ = Ea, we have 

1 /■ ,„ 1 



— chdE—f dE' AgaiE',E,K\K) 
e2|e|Z''g°^(a|a'i?''=(0)|n)(n|a'|a 



2M ^ \ Ea-En 

kl /fy-/ 



{a\a'\n){n\aW^HO)W)\r^ku^, 



^ N I I ^ y_JL^ ) [Bf (K' , K) + ^2 (K' , K) 



M 2^J_^ ^ \{Ea~En){ea~UJ-Enil-iO)) 



{a\a^\n){n\aW''{uj)\a) \(jki^j^f ^g^^ 



(Ea - £n)(ea + W - " iO)) 

where 

C"'(K',K) = Sf(K',K)-B2"'(K',K) = -^efe,,7^^<5(K'-K) (86) 

and ekis is the Levi-Civita symbol (ei23 = £312 = £231 — 1, £321 = £132 — £213 = ~1, and £^(5 = if at least two of the 
indices are equal each other). It can easily be shown that the term containing [Sj^' + i?!'] in equation (|8^) is equal to 
zero. Indeed, using the identity 

a'^i[Hu,x'], (87) 

where Hd = —ia ■ V + f3m + V , we obtain 

{a\a'D'''{0)\n){n\a^\a) {a\a^ \n) {n\a' D''' {0)\a) \ 



= i J2 i-{a\aW'{0)\n)(n\x^\a) + {a\x^\n){n\a'D'''{0)\a)) 

n 

= i{-{a\aW''{0)x^\a) + {a\x^ aW' {0)\a)) = . 

Here we have taken into account that for the case of an extended-charge nucleus considered here there are no degenerate 
states of different parity. Taking into account that to zeroth order the denominator in equation ( p3| ) is equal to 1 , we 
have 

e^\e\Zi ^^i f°° ^ '"J^'i {a\aW''{Lu)\n){n\a'\a) 



(89) 



-{Sa - £n){£a - UJ - En{l ~ iO)) 

{a\a'\n)(n\aW'''{uj)\a) 



[ea - £n){£a + - e«(l - iO)) 
Using identity (p7|), we obtain 
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/oo 
dcu {{a\x'G{uj)D''{ea - ij)\a) + {a\D''{ea - u;)G{oj)x'\a)} , (90) 
-oo 

where D^{u!) is defined by equation ([n]). 

The other contributions to AiJ^^^'^defined by equation (^6|) are derived in the same way. As a result of this 
derivation, keeping only the state-dependent contributions, we have 

/oo 
dio {(a|a;'GH(ma'= -/)|a) + {a\{ma'' - p^)G{u)x'\a)} . (91) 
-oo 

The sum of expressions (1^) and d^) gives the correction /\E^'^^\ For the correction AiJ^^) = A£;(2'») + AE('^^\ 
where AE^'^"''' is defined by equation we obtain 



Ai?(2) = / duj {{a\x'G{cu + £„)[/ - ^?'(^)]|a> 

2M 2tt 

+ (a|[/ - D''iLo)]Giu + ea)x'\a)} . (92) 

For practical calculations, it is convenient to represent this correction by the sum of the lower-order and higher-order 
terms, 

AE^^^ = AE^'^ + AE^'^ , (93) 
AE^^^ = -^ekisnma\x'[p' - i?'=(0)]|a) , (94) 

/OO 1 

-(a|[/ - D''iu;)]G{LO + Sa)a'\a)} . (95) 
Here, as in equations (|76|)-(|7^), the term AE^j^'^ contains the non-relativistic and lowest-order relativistic contributions 

(2) 

and a part of the higher-order relativistic contributions while the term AE^^ contains all the higher-order relativistic 

(2) 

contributions which are not included into AE^^ . 

The total energy shift of first order in m/M and of first order in the interaction with the classical magnetic field is 
^^(tot) ^ A£;(i) -I- A£;(^) . In the compact form, AE^^'^ and A£;(2) are defined by equations (|7|)and (||) , respectivelv. 
For practical calculations, it is more convenient to use the representation given by equations (|76|)-(|7q) and (|93|)-(P5[). 
To derive these equations, we have assumed that V{'x.) deviates from the pure Coulomb potential. However, since the 
final formulas for the energy shift exhibit neither singularities nor ambiguities when V(x) — > Vc = — aZ/|x|, the pure 
Coulomb potential can be used in these formulas as well. 

The corresponding correction to the bound-electron g factor is defined as 

A£'(*°*) 

Ag=——, (96) 

where fiQ = |e|/(2m) is the Bohr magneton and rrij is the angular momentum projection of the state under consider- 
ation. Here and below we assume that Ti. is directed along the z axis. 

V. RECOIL CORRECTION TO THE BOUND-ELECTRON g FACTOR TO LOWEST ORDERS IN 

aZ 

To the lowest-order relativistic approximation, the recoil correction to the g factor is given by 

where AE'^^-' and AE^'' are defined by equations ( [tt] ) and (^4|), respectively. Let us calculate this correction for an 
arbitrary state of a hydrogenlike atom. For the case of the point-charge nucleus that we will consider, this calculation 
can be performed analytically. 

Consider first the calculation of AE^\ According to equation (|77|), it is 



Ae['' = i^^-^^l [p' - • P + (« • n)(n . p) \a) , (98) 
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where n = r/r. Taking into account that = (a • p)^ and (a • p) = 7?d — Pm — Vc, one easily obtains 

(<5a|p2|a) = {Sa\{ea + /3m - yc)(£a - /3m - yc)|a) 
+i((5a|(a- Vyc)|a). 

The second term in equation ( p8| ) can be transformed as (see, e.g., ||35[| ) 

- (Sal [a • p + (a • n)(n • p)]|a) = -{Sal 2ea - 2(3m - 2Vc 

r r I 

+ • n)(/3K + 1) \a) , 



(99) 



(100) 



where k = (^l)"'^'^^''^(i + 1/2) is the relativistic angular quantum number of the state a, j is the total angular 
momentum, and I — jztl/2 defines the parity of the state. The wave function correction \Sa) defined by equation ( [75| ) 
can easily be found analytically using the method of the generalized virial relations for the Dirac equation developed 
in p6[ |. Since the operator sandwiched between \a) and \Sa) in equation for AE^^"^ conserves the angular quantum 
numbers, we need only that component of \Sa) which has the same angular quantum numbers as the unperturbed 
state I a). This component is 



\Sa) 



X(r)f]„„^-(n) 
iY{r)n_^rnj (n) 



where 



Xir) = bo 



2m,K — m + 2Kea aZ 



Y{r) 
bo 



= br 



2m? m? 
2mK + m — 2Kea aZ 



2rr? 



fir) 



2m2 
2m2 



air) , 

fir), 



2-^j(j- + l)-'^' 

g{r) and /(r) are the radial parts of the unperturbed wave function defined as 

g{r)Vt^yrL, (n) 
if{r)n^i^,n, (n) 



Integrating over the angular variables in equations (M) and (llOOf), we find 



^Ei'^ ^ 17 y„ r2|x(r)g(r)[-2Fcm - + el - m"] 



+Y{r)f{T)[2Vcm - + el - 
+ [X{r)f{r)+Y{r)g{r)]^K}. 



Substituting expressions ( |102| ) and ( |103| ) into equation (|106| ), we obtain 

j^^ 2A.£,-m ^o + (c,Z)2-C-i + [el - rr?)-C^ 
L m m m 

-aZ- 



^ M 



-{aZf 



where we have used the notations 



2m? 



A- 



i^a-^')^,B' + ocZ 



3m,K — 2K?ea 
2m^ 



B 



A" = 



dr r'+^ig' + f'). 



B" = 



dr 



C ^2 



dr 



gf- 



(101) 

(102) 

(103) 
(104) 

(105) 



(106) 



(107) 



(108) 



The integrals A^ , B'^ , and C are easily calculated by the recurrent equations given in (the relevant equations 
were first derived in p^). We have 
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1 _ [aZYnm 



9 9 

re m — et 



2(aZ)3[2K(7 + ?v)-7V] 



^1 _ 2K£a - TO 



2™^ ' 



7V4(472 - 1)7 



5° = ^ 

TO 



7V37 
TO^ — 
aZm ' 



where is the radial quantum number, TV = \/ n'^ — 2nr(|K| — 7), and n = n,. + |k| is the principal 
Substituting formulas (109) into equation (107), we obtain 



(1) 



2 M 2j{j + 1) to2 



Consider now the term AE'P-' . According to equation ([Ml) , we have 



(2) 



e n 



■L - 2 X p]. - [r X D(0)],)|a) 



Integrating over the angular and radial variables, we obtain 



aZ , 



L o ;i,f o„v,-_|.i)\-^^-^ ^ ^ 4 ^ ^ 2to 



2 M 2j(j + 1) 
+ (2j-0(2i-/ + l) 



2m 



, TO — e 



For the sum AE'l = A£'[^' + AEl"' , we find 



p(2) 



The related contribution to the g factor is 
To the two lowest orders in aZ , we have 



TO 1 

A5L = -T7 



. ^ e H 2k^£^ + KTOEa — TO^ 

Ai^L = ^ rn-i 

2M 2m2j(j + l) ' 

m 2k^£^ + KTOffa — m? 



M 2m'^j{j + 1) 



M jU + 1) 



2 2 V 4/ 



2 , ("^)^ 



For an ns state, formula (114) yields 



TO 2 2e^ — mSa — to^ 
M 3 ^ 



while formula (|11E| ) gives 



AffL 



TO (aZ)2 
M n2 



Formula (117) agrees with the related result obtained for the Is state in ||20|-|2 
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VI. CONCLUSION 



In this paper we formulated the systematic QED method for calculations ol the nuclear recoil corrections to the 
energy levels and to the electronic g factor in atoms to first order in m/M and to all orders in aZ. Employing this 
method, we derived the complete aZ-dependence expression for the recoil correction to the bound-electron g factor 
in a hydrogenlike atom (equations (|73|), B2[)) . This expression was also represented as the sum of the lower-order and 
higher-order terms (equations ([76|)-(]78|), (|93|)-(95)). The lower-order term contains all the non-relativistic and lowest- 
order relativistic contributions and a part of the higher-order corrections. For an arbitrary state of a hydrogenlike 
atom, a simple analytical formula was derived for this term (equation ( |ll3| )). As to the higher-order term, we expect 
that its numerical evaluation can be performed in the same way as the evaluation of the related contribution to the 
energy levels [|8|-|0). This calculat ion, which is equally important for low- and high-Z systems, is under way and will 
be pubhshed elsewhere. 
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